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^ ■ Abstract 



We study generic single- field dark energy models, by a parametrization of the most general theory 



\ of their perturbations around a given background, including higher derivative terms. In appropriate 

CN . limits this approach reproduces standard quintessence, fc-essence and ghost condensation. We find 
00 ■ 

■ no general pathology associated to an equation of state wq < —1 or in crossing the phantom divide 
wq = — 1. Stability requires that the wq < — 1 side of dark energy behaves, on cosmological scales, 

00 \ as a /c-essence fluid with a virtually zero speed of sound. This implies that one should set the speed 

' of sound to zero when comparing with data models with wq < — 1 or crossing the phantom divide, 

k J ' summarize the theoretical and stability constraints on the quintessential plane {1-\-wq) vs. speed 

ij ■ of sound squared. 

d, 

1 Introduction 

The origin of the present acceleration of the Universe is likely to be the most important theoretical 
problem in physics today. Given the general reluctance in accepting as explanation an incredibly 
small cosmological constant and the absence of compelling alternatives, it seems that one should keep 
an open-minded approach, concentrating on very general theoretical constraints and on observables 
more than on specific models. 

In this paper we study in generality, and focusing on perturbations, dark energy scenarios where 
the dark sector is described by a single scalar degree of freedom, without direct coupling to matter 
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(in the Einstein frame). We will often call this general model quintessence, although in the literature 
this name is usually reserved to a scalar field with a canonical kinetic term. 

Following [1, 2], we will rewrite the scalar field Lagrangian in order to make explicit what is the 
most general theory of quintessence perturbations around a given background solution characterized 
by its pressure pq and energy density pq. The freedom that we have after the unperturbed history 
is fixed is made clear in this way. This separation is particularly important given that a host of 
new experiments is going to test dark energy clustering properties [3] . In our formalism, the general 
theoretical constraints on single field models are also made clear. 

In particular, we will study whether a single field model that is safe from ghost and gradient 
instabilities can have an equation of state wq < —1, where wq = pq/pq- In this regime, the 
stability of the model can be guaranteed by the presence of higher derivative operators, a conclusion 
already reached in [1], where single field models were studied focusing on the constraints enforced 
by stability. Here, after reviewing and extending the results of [1], we will concentrate on the 
behaviour of cosmological perturbations, which are relevant for observations. On cosmological scales 
we find that these higher derivative terms are irrelevant for the phenomenology, so that a model with 
wq < — 1 simply behaves as a A;-essence fluid with virtually zero speed of sound. Higher derivative 
terms are relevant for cosmology only when the equation of state gets very close (and experimentally 
indistinguishable from) a cosmological constant. In this limit our general Lagrangian reduces to 
the Ghost Condensate theory [4] smoothly connecting quintessence to this theory of modification of 
gravity. Notice that, as detailed in Appendix A, we are interested in a regime where higher derivative 
terms do not introduce additional degrees of freedom (contrary to what happens for example in [5]). 
We find it convenient to summarize our results in the plane (1 + wq)Qq vs. c^, where $7q is the 
quintessence contribution to the critical density. We dub this plane of parameters the quintessential 
plane. 

We also study the issue of whether it is possible to cross the so-called phantom divide wq = — 1 
[6, 7]. We find that the speed of sound vanishes exactly at the divide [8, 9] and since quintessence 
may remain stable for wq < — 1 there is no general pathology associated with the crossing. We show 
this explicitly with an example. The phantom divide can be crossed with a single scalar degree of 
freedom, without introducing ghost-like fields. 

The paper is organized as follows. In section 2 we study the most general theory of single 
field quintessence, taking into account higher derivative operators and focusing on the stability 
constraints following [1]. In section 3 we study the phenomenology in various limits, considering also 
the gravitational efi'ect of dark matter on quintessence. In section 4 we consider the issue of crossing 
the phantom divide wq = — 1 and we show explicit examples of the crossing without pathologies. 
In section 5 we concentrate on another kind of higher derivative operators [1, 2], different from the 
ones studied for the Ghost Condensate. Although the phenomenology on cosmological scales does 
not change, the modification of gravity at short distances is quite different. Conclusions are drawn 
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in section 6. 

Several issues concerning our effective theory approach are left to the appendix. Appendix A is 
devoted to reviewing how higher derivative operators must be treated in the effective field theory 
approach. In appendix B we write down the stress-energy tensor for the action discussed in the main 
text. In appendix C we derive the most general action for quintessence perturbations following the 
approach of [1, 2]. Finally, in appendix D we discuss the modification of gravity induced by the kind 
of higher derivative operators that were not studied in [4] . 

2 Effective theory of quintessence 

Our aim is to study the most general theory of quintessence perturbations. We will do it step by 
step, first by considering a model with an action containing at most a single derivative acting on 
the field. This is known as A;-essence [10, 11] and it will be possible to write the action for the 
perturbations in such a way as to make explicit the dependence on the background energy density 
and pressure pq and pq. Then we will add higher derivative operators to the A;-essence action in such 
a way as to leave the background invariant. In this section we will consider the kind of operators 
introduced in the context of ghost condensation. Other higher derivative operators will be discussed 
later in section 5. 

An alternative derivation of the most general action for quintessence perturbations is given in 
appendix C following the approach of refs. [1,2], that consists in writing down all the terms preserving 
the symmetries of the system in a 'unitary' gauge, where the quintessence perturbation is set to zero 
and appears as a scalar metric degree of freedom. This approach is elegant and straightforward but 
less pedagogical than the one adopted here. Both approaches lead to the same physical results. 

2.1 The limit of /s-essence 

Let us start with a /c-essence action 

S =y"d4 xV=5 P{<l>. X), X = -g>'^d^<j)d,<i> . (2.1) 

We assume a flat Fricdmann-Lemaitre-Robertson- Walker (FLRW) Universe with metric ds^ = 
-dt^ + a?{t)<lx^. Initially, we will treat this as a fixed background and neglect the perturbations of 
the metric. 

To describe perturbations around a given background solution 4>o{t), it is useful to write the 
scalar field as 

= 0o(i + 7r(i,f)) , (2.2) 

and expand the action (2.1) in terms of tt. In the following, we are going to assume that the function 
0o(i) is strictly monotonic, 0o(O 7^ 0, to avoid singularities in the relation between ^ and tt. 
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Using the expansions 

<p{t,x) = 

X{t,x) = Xo + XoTT + ^XoTT^ + 2Xo7r + 2Xo7r^ + Xo ( - 



1 



(Vtt) 



+ ... , 



(2.3) 
(2.4) 



where Xq = (pQ , we have, up to second order in tt, 

1 ^ 



S 



Po + Pott + ^tPott^ + 2PxXo7r + 2 {PxXoYTrTr + PxXo - 



(Vtt)^ 



+ 2PxxX^7r^ 
(2.5) 

where Px = dP/dX\o and Pxx = 5^P/5X^|o. The term Pq does not affect perturbations as it is 
independent of tt, while one can verify that the linear terms cancel using the background equation 
of motion. Indeed, by integrating by parts the term tttt and making use of the background equation 
of motion, after some manipulations we are left with 



■I 



d^xa^ 



{PxXo + 2PxxXi) TT^ - PxXo^^ + SHPxXo tt^ 



(2.6) 



We can now rewrite the coefficients of this expansion in terms of the stress-energy tensor of the 
background solution. Prom the definition of the stress-energy tensor, 

2 6S 



one obtains the background energy density and pressure, 

PQ = 2XoPx - Po , PQ=Po- 
Using these expressions the action above can be cast in the form 

S=jd^xa^ ^ [pQ +PQ + AM^) TT^ - 1 {pQ + pq) + ^H{pQ + PQ^ 



(2.7) 



(2i 



(2.9) 



Here we have defined = PxxX^^ where M has the dimension of a mass. 

At this stage wc can straightforwardly introduce the coupling with metric perturbations. This 
coupling is particularly simple in synchronous gauge where the metric takes the form 



ds^ = -dt^ + a^{t){6ij + hij)dx^dxP . 



(2.10) 



Indeed, at quadratic order in the action the coupling with gravity only comes through the perturbed 
-y/^ in the action. Replacing with a^(l -|- h/2) in eq. (2.5) we have 



:/d^.a3(l4 



Pqtt + 2PxXo7: + 



(2.11) 



4 



Integrating by parts and using again the background equation of motion one gets the full action for 



(2.12) 



The quadratic Lagrangian is thus specified by the functions (/9q and M^(t). It is 

important to stress that these two functions are completely unconstrained. For any choice of these 
two functions one can in fact construct a Lagrangian P{<p,X), such that the quadratic Lagrangian 
around the unperturbed solution has the form (2.12). Let us see explicitly how this works. First 
of all, given pg + pg (and possibly the other contributions to the total energy density and pressure 
coming from other components) one can find the two functions PQ{t) andpQ(i) solving the Friedmann 
equations.^ At this point it is easy to check that the Lagrangian 

Picf>,X) = ^{pQ - pq)(0) + ^(pq+pq){cP)X + ^M\cP){X - If (2.13) 

has the solution <p = t, gives the requested pressure and energy density as a function of time, and 
gives eq. (2.12) as the quadratic action for perturbations. Note that the dimension of (f) is that of 
an inverse of a mass. The coefHcient is time dependent and for quintessence we expect that it 
varies with a time scale of order H~^. Somewhere in this paper, to simplify the calculations we take 
= const. 

One advantage of the action (2.12) is that the coefficients of all terms are physically measurable 
quantities - we will see below that is related to the sound speed of perturbations. This standard 
form of rewriting the action of quintessence perturbations does not suffer from field redefinition 
ambiguities. Indeed, there is an infinite number of physically equivalent Lagrangians P{4>, X) related 
by field redefinitions (['{4') but they all give the same action (2.12). Note that for most purposes 
the explicit construction of the action in terms of cp is irrelevant. Indeed, one is free to choose any 
function pgia) to describe the evolution of the quintessence energy density as a function of the scale 
factor. Then the action (2.12) will describe the perturbations around this background. In particular, 
one can always make a field redefinition such that = t, as we did in eq. (2.13). 

Let us see what are the theoretical constraints that we can put on the general form of the action 
(2.12). A basic requirement that we will impose on our theory is that it is not plagued by ghosts, 
i.e. that its kinetic energy term is positive, 

PQ+PQ + 4M^>0. (2.14) 



^One has to solve the continuity equation for quintessence, PQ + iH{pQ-\-pQ) = 0, where Pq+Pq is a known 
function of time, together with the Friedmann equation, H'^ = (pg + /Orest)/3M|[ where Prest includes all the 
additional sources of energy density in the Universe. These equations can be integrated up to a constant in 
the initial condition. This ambiguity corresponds to a shift in the cosmological constant which does not enter 
in the action for perturbations. 
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The presence of a sector with the "wrong" sign of the energy impUes that the Hamiltonian is 
unbounded from below. If one studies this sector alone, no pathology arises as the sign of the energy 
is a matter of convention. However, quintessence is (at least) gravitationally coupled to the rest of 
the world, so that there is the danger of exchanging energy between a healthy sector and a negative- 
energy one without bound. Classically this is not a problem if quintessence perturbations are very 
small and remain in the linear regime. Therefore, for a quintessence with negligible clustering (with 
speed of sound Cg ~ 1), there is no obvious classical danger. At the quantum level the situation 
is more pathological. The vacuum is unstable to the spontaneous decay into positive and negative 
energy states and the decay rate is UV divergent and strictly infinite in Lorentz invariant theories [12]. 
Although it has been shown that it is possible to cut-off this instability in a non-Lorentz invariant 
theory [13], in this paper we take a more conservative approach and forbid the existence of ghosts. 

If we set = in the general action (2.12) we reduce to the case of a standard quintessence 
field with canonical kinetic term (50)^. In this case, forbidding the ghost implies PQ+PQ > 0. Thus, 
as it is well known, a scalar field with minimal kinetic term can violate the null energy condition, 
equivalent in a cosmological setting to p + p < 0, only if it is a ghost [14]. In this simple case the 
speed of sound of scalar fluctuations is = 1. When does not vanish the speed of sound of 
fluctuations differs from unity [15] and reads 

2 _ PQ+PQ .2 15^ 

'~ PQ+PQ + AM^ ■ ^ ' 

One can see that, for PQ + pq > 0, i.e. for positive c^, < implies that scalar perturbations 
propagate super-luminally [16, 17]. This is problematic in a theory with a Lorentz invariant UV 
completion [18]. 

Prom the action (2.12) we see that in the presence of there is no generic connection between 
the violation of the null energy condition and a wrong sign of the time-kinetic term. The coefflcient 
in front of tt^ can be positive also when pq+pq < 0. On the other hand, pq + pq fixes the sign of 
the term in front of the spatial kinetic term (Vtt)^ [19]. Thus, in absence of ghosts, the violation of 
the null energy condition implies a negative speed of sound squared. The constraints that we have 
derived are summarized in the quintessential plane 1 + wq vs. c^, represented in figure 1. 

An imaginary speed of sound (c^ < 0) represents a gradient instability of the system. Taking 
^ \PQ + Pq\ gradient term is suppressed and the instability is irrelevant for scales of 
cosmological interest in the sense that the instability time is much longer than the age of the 
Universe. Still, the instability is relevant for short wavelengths so that it seems difficult to make 
sense of the wq < — 1 region [8]. However, this conclusion is too hasty. Indeed, when the term 
(Vtt)^ is suppressed, higher derivative operators may become relevant as we will discuss in the next 
section. 
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Figure 1: The quintessential plane 1 + wq vs. Cg in the case of k-essence. If we require the absence 
of ghosts, the sign of the spatial kinetic term is fixed to be the same as 1 + wq, so that one has to 
worry about gradient instabilities for 1 + wq < 0. For 1 + wq > one has superluminal propagation 
for < 0. 

2.2 Higher derivative terms and stability 

In the previous section we saw that in the limit in which the quintessence speed of sound goes to 
zero, the standard spatial kinetic term vanishes. Obviously the action (2.1) is not the end of the 
story: the full Lagrangian will contain higher derivative operators such as Q{X)R{\I\(j)) and these will 
give rise to the leading higher derivative spatial kinetic term. A generic higher derivative operator 
will however change the background solution of (2.1), while here we want to study the effect of 
higher derivative operators on quintessence perturbations around a given background, as we did in 
eq. (2.13). To keep the background unchanged, let us add to the Lagrangian (2.13) the operator 

CM = - — {a(t> + 3H{(t>)f. (2.16) 
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For reasons that will become clear later, we need > 0? This term does not alter the background 



evolution 4> = t, PQ{t) and PQ{t). Indeed, □(?!) + 3ii"(0) vanishes on the background so that the oper- 
ator is explicitly quadratic in the perturbations. At quadratic order this operator reads (neglecting 
for the moment metric perturbations) 



One may worry about the presence of terms with higher time derivatives, as these would naively 
be associated with additional solutions of the equation of motion. However, if one compares M^tt^ 
with the standard time kinetic term M^tt^ of eq. (2.12) - assuming M ~ M - the former is always 
suppressed with respect to the latter for frequencies below the scale M. In general, wc expect 
that for frequencies a; ~ M all operators containing higher time derivatives become important, so 
that the scale M ~ M sets the maximum energy scale for which the theory makes sense: it is the 
energy cutoff. This is the standard situation in an effective field theory: higher derivative terms 
become important for energies of the order of the cutoff and at lower energies they must be treated 
perturbatively. In particular, there is no physical meaning in the new solutions that arise from 
taking higher and higher time derivatives. We postpone a complete discussion about this point to 
appendix A. Notice that the same argument cannot be used for the operator — M^(V^7r)^/2. Indeed, 
in the limit of small pq +PQ there is no spatial kinetic term of the form (Vtt)^ so that — M^(V^7r)^/2 
is the leading spatial kinetic term. At short scales we have a non-relativistic dispersion relation of 
the form lo ~ k'^/M which implies that energy and momentum behave very differently (as we will 
see in appendix A they have different scaling dimensions). In particular, when comparing the first 
and last terms in the brackets in eq. (2.17) we have V^tt ~ Mtt S> tt, for energies below the cutoff. 
This means that we can drop it altogether from eq. (2.17). 

It is important to stress that there is no fine tuning in the limit \pQ +pq\ <^ - or equivalently 
\(?g\ <^ 1 - as this limit is "technically natural", i.e., there is a symmetry that is recovered in the 
limit PQ + PQ = 0. Indeed, as shown below, in this limit we obtain the Ghost Condensate theory 
[4], which is invariant under the shift symmetry ^ — + A. In the presence of this symmetry the 
expansion of the Universe drives the background solution to pQ + pQ = 0. Thus, models with very 
small speed of sound should be thought of as small deformations of the Ghost Condensate limit 



Let us come back to the issue of stability for pg + pg < including the new higher gradient 
term (2.17) and follow the discussion in [1]. As we discussed in the previous section, the only 
dangerous modes are those on scales much smaller than the Hubble radius, as their instability rate 
can be arbitrarily large; wc thus concentrate on k/a^ H. In this regime the operator (2.17) further 
simplifies as the second and third terms are negligible with respect to V^tt. Considering the action 

^As also can have a time dependence on a time scale of order H~^. For simplicity, in the following 
we assume that it is constant. 




(2.17) 



[20, 1, 2]. 
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(2.12) with the addition of the only remaining operator — M^(V^7r)^/2, the dispersion relation of tt 
is thus modified to 2 4 

[PQ +PQ + 4M4) u;^-{pq+pq)^-M^'^ = 0. (2.18) 

For Pq +Pq < the A;^ gradient term has the unstable sign, but in the presence of the new operator 
this instability is confined to sufficiently large scales. In particular the fastest rate of instability is 
given by 

^Ld--^^^^^, (2.19) 
where we have taken 3> \PQ + Pq\- This gradient instability is not dangerous when it is slower 
than the Hubble rate, i.e., when 

_PQ+M<H (2.20) 

It is clear that a larger M makes the gradient instability slower. 

However, a large M sources another form of instability, which contrarily to the gradient instability 
is already present for PQ + PQ = and was originally discussed for the Ghost Condensate theory 
in [4]. Indeed, when the coupling with gravity is taken into account, the system shows a sort of 
Jeans instability similarly to a standard fluid. To see this, let us introduce metric perturbations 
and consider the limit PQ + PQ = 0. In this case, the complete Lagrangian reads, neglecting the 
expansion of the Universe, 

2- 

(2.21) 



S = ld*x 



2M^7r2 




We see that a large M enhances the mixing of tt with gravity, i.e. the Jeans instability. 
Thus, the equation of motion for tt reads, in this case. 

The gravitational perturbation h is sourced by the perturbations of tt through Einstein equations. 
In synchronous gauge, h satisfies [21] 

h = -^iSPQ + ^SPQ), (2.23) 

where we have neglected the expansion of the Universe, and we have introduced the reduced Planck 
mass Mpj = (SttG)"-*^. The stress-energy tensor can be straightforwardly derived by varying the 
action (the complete expression of the stress-energy tensor is given in appendix B) and the leading 
term is Spg = AM^-k while the pressure perturbation is negligible.^ The solution of eq. (2.23) can 



^Indeed, we have 




V^ttJ , (2.24) 

(2.25) 



be plugged back into eq. (2.22). This yields the equation of motion of tt taking into account its 
gravitational back-reaction, 

7f + — — VV = V\. 2.26 

Mixing with gravity induces an unstable fc^ term in the dispersion relation, similarly to the gradient 
instability discussed above. We can compute again the fastest instability rate. 



.2 

Jeans 



Ml, 



(2.27) 



As expected, in this case the instability gets worse for large M^, i.e., when the mixing with gravity 
is enhanced. By imposing that this instability rate is smaller than the Hubble rate^ we obtain 

MM2 



Ml, 



<H. (2.28) 



Requiring that both stability conditions (2.20) and (2.28) are satisfied we get the window [1] 

MM2 

-(l + WQ)nQ<j^^<l. (2.29) 

In conclusion, considering higher derivative terms, a quintessence model with wq < — 1 can be 
completely stable inside the window of parameters (2.29). On the other hand, eq. (2.29) indicates 
that it is difficult to avoid instabilities when (1 + wq)CIq ^ — 1. These stability constraints were 
already obtained, for J7q = 1, in [1]. 



3 Phenomenology on the quintessential plane 

3.1 /c-essence vs. Ghost Condensate 

Coming back to the quintessential plane of figure 1, in the previous section we have learned an 
important lesson: the gradient instabilities for wq < — 1 can be made harmless by higher derivative 
operators. Thus, part of the lower left quadrant of the quintessential plane is allowed. 

To discuss the phenomenology of these models (for a related discussion see [23] ) let us write the 
full action for perturbations including the higher derivative operator (2.16): 



-I 




2- 

(3.1) 



2 a 



2 



For k <C M, M^V^n <^ M^tt. Moreover, eq. (2.23) shows that also the h terms can be neglected in front of 

M^TT, SO that the operator proportional to gives a negligible contribution to the stress-energy tensor. 

more careful analysis [22] indicates that this condition is very conservative and much larger instability 

rates can be experimentally allowed. 
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First of all, note that it is not possible to switch off quintessence perturbations for PQ + PQ 7^ 0; 
doing it by hand would give gauge dependent unphysical results. On the other hand, the converse is 
not true: even for Pq+Pq = ^ perturbations may still be present, as in the Ghost Condensate case. 

We saw that the operator in the second line of eq. (3.1) allows the stabilization of the short 
scale gradient instability; on the other hand, for cosmological purposes we are interested in very 
large scales. Let us see whether this operator is relevant for scales of the order of the Hubble radius 
(although our discussion will extend to all scales of cosmological interest). We want to show that 
when 

|/9q+pq|»m2^2 (3.2) 

the higher derivative operator can be neglected when discussing the cosmological clustering of 
quintessence. In this case we reduce to a standard fc-essence model, with the only difference that 
there are no short-scale instabilities even for wq < — 1. On the other hand in the opposite case, 

\pq+Pq\^M^H\ (3.3) 

all the terms in the action (3.1) proportional to pq + pq can be neglected. In this case the model 
reduces to the Ghost Condensate theory. 

Verifying the existence of these two regimes is quite straightforward. For instance the dispersion 
relation at k/a ^ H \s dominated cither by {pq +pQ)(V7r)^ or by M^(V^7r)^ depending on the 
hierarchy between \pq + pq\ and M^H^. The same applies for the operators involving the metric 
perturbation, {pq + pq)h'K and M'^hS/^'K. This check can be done for all the other operators, by 
taking V/a H and considering that time derivatives axe at most of order H. The existence of 
these two regimes can also be seen by looking at the stress-energy tensor (see appendix B). 

Now we can go back and complete our quintessential plane. When wq is close to —1, 

-^S(l+w«)ng<^. (3.4) 

the model behaves as the Ghost Condensate. We can estimate the width of this region by imposing 
the absence of Jeans instability, eq. (2.28). Assuming M ~ M one gets a rough upper bound: 
M < 10 MeV [4]. A more accurate analysis shows that this limit is much too conservative and it 
can be relaxed to M < 100 GeV [22]. Even in this case the window above is extremely tiny: 

\1 + wq\nq< 10-^'' . (3.5) 

We can therefore draw an important conclusion: only for values of wq which are observationally 
indistinguishable from the cosmological constant, does quintessence behave as the Ghost Condensate 
on cosmological scales. This regime corresponds to the strip around the horizontal axis wq = — 1 in 
figure 2. Notice that in this region the dispersion relation is of the form uj oc k"^, so that the speed 
of sound is not well defined, i.e. it becomes scale dependent. 
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(1 + w;q)1)q 
A 




Figure 2: On the quintessential plane, we show the theoretical constraints on the equation of state 
and speed of sound of quintessence, in the presence of the operator M . Instability regions are dashed. 
Where 1 + wq and have opposite sign we have a ghost-like instability corresponding to negative 
kinetic energy. For wq < —1, the dashed regions in the left-lower panel is unstable by gradient 
(c^ < —HM/M'^) and Jeans ((1 + wq)^q < —1) instabilities, while the strip close to the vertical 
axis corresponds to the stability window (2.29). Furthermore, the strip around the horizontal axis 
given in eq. (3.4) corresponds to the Ghost Condensate. Above this region, in the right-upper panel, 
we find standard k-essence. 

On the other hand, for any value of wq which is appreciably different from the one of the 
cosmological constant, the model reduces to fc-essence as higher derivative terms are cosmologically 
irrelevant. Their only role is to stabilize the short scale gradient instabilities for wq < — 1. Although 
in practice not relevant, note however that wq cannot be made arbitrarily negative. This is shown 
by eq. (2.29) and, in the quintessential plane, it excludes the bottom shaded region of the lower-left 
quadrant. 

Let us now constrain the values of the speed of sound c^. For wq > — 1 there are no constraints, 
besides the possible limit < 1 already discussed. For wq < — 1 the speed of sound is negative and 
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very small, as it is constrained by the absence of gradient instability, eq. (2.20), 

2^, PQ+PQ ^HM 

We can numerically constrain the right-hand side of this equation by considering that the scales 
M ^ M represent the cutoff of our effective field theory. By requiring this cutoff to be larger than 
the minimum scale at which gravity has been proved, i.e., M > 10~^eV, and using in eq. (3.6) the 
value of the Hubble parameter today, Hq ~ (10~^eV)^/Mpi, we obtain 

n \ V2 

T7^) ~10-=^°. (3.7) 

Thus, for all practical purposes the speed of sound can be taken to be exactly zero. On the quintessen- 
tial plane in figure 2, in the lower- left quadrant, we can only live in a tiny strip along the vertical 
axis. Notice however that there is no fine tuning in keeping extremely small. Indeed, as we 
discussed, in the limit of Ghost Condensate vanishes exactly for symmetry reasons. Thus, the 
speed of sound remains small for tiny deviations from this limit. 

3.2 Including dark matter 

After the discussion about the stability constraints, we would like to understand the dynamics of 
quintessence perturbations and their impact on cosmological observations. In order to do this, 
we will now study quintessence in the presence of cold dark matter, which gravitationally sources 
quintessence perturbations. A thorough analysis of the phenomenology of these models is beyond 
the scope of this paper. Here we want to focus on the main qualitative features in the various limits. 

Let us start from the Ghost Condensate limit (3.3). It is known that the Ghost Condensate 
affects only short scales, i.e., tt perturbations induce a modification of the Newtonian potential at 
scales parametrically smaller than the Hubble scale [4]. Therefore, we expect to have extremely 
small effects on cosmological scales. To verify that this is the case, we can study the action (3.1) in 
the limit of Pq+Pq = 0. This reads 



2MV-^ 3ifi-3i7. + ^-^ 
2 \ 2 

For simplicity, let us momentarily disregard the first two terms in parentheses, 



(3.8) 



S 



2M^7r2 - 



f h V^ttV 



2 \ 2 a 



2 



(3.9) 



Notice that this is the action used in the Ghost Condensate paper [4]. The equation of motion for 
the TT perturbations is given by 

13 



The gradient term on the left hand side can be neglected on cosmological scales. Indeed, the time 
derivatives will be at least of order ~ MA;^/(aM)^, so that the friction term in the previous equation 
will always dominate the gradient term fov k/a ^ H. As we want to show that Ghost Condensate 
perturbations remain small, we assume that the dark matter dominates the perturbed Einstein 
equations. The validity of this assumption can be checked a posteriori. 

In a matter dominated Universe h = —25m [21], with 5^ = ?>Pml Pm-, which, using = HSm 
and the background Priedmann equation, leads to 

A = -A^. (3.11) 

We can now replace this as the source of Ghost Condensate perturbations on the right-hand side of 
eq. (3.10). This yields, neglecting the gradient term, 

TT + 3il TT 2i ^ ttt;— ■ (3.12) 

12M4m2j a^H ^ ' 

If we now assume that the initial quintessence perturbations arc small so that the homogeneous 
solutions arc sub-dominant, similarly to what happens in standard quintessence [24], this equation 
can be solved to give 

24M4Af|j a2if2 • 

As we discussed in the previous section, the energy density and pressure perturbations of the 
Ghost Condensate are dominated by the operator so that bpq ~ AM'^ir and Spg 2± 0. Thus, on 
cosmological scales, 

^PQ ~ TTT^Pm ■ (3-14) 

Prom the simple estimate of M below eq. (3.4), we conclude that quintessence perturbations are 
negligibly small with respect to dark matter perturbations, 6pQ < 10~^'^5pm- It is straightforward 
to generalize this analysis including the two terms in parentheses of eq. (3.8) previously neglected 
and verify that eq. (3.14) remains valid. The conclusion of eq. (3.14) is quantitatively consistent 
with the (small) modification of the Newton law derived in [4] , as one can check for example in their 
eq. (7.11) for k/a ~ H and uo ^ H. 

Close to the wq = — 1 line, we saw that there are no appreciable effects of perturbations on 
cosmological scales which can help in distinguishing quintessence from a cosmological constant; all 
the interesting dynamics is limited to short scales. As we move away from the wq = — 1 line we 
enter in the A;-essence regime, as we pointed out in the previous section. The case wq > — 1 is well 
studied in the literature (see, for instance, [25, 26]). The case wq < — 1 is much less studied: here 
we have a negative speed of sound squared that is so small - see eq. (3.7) - that can be taken to be 
zero for all practical purposes. With a small speed of sound we expect quintessence to cluster on 
scales shorter than the Hubble radius driven by dark matter gravitational potential wells. To study 
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this, let us repeat the calculation we just did in the Ghost Condensate case for a fc-essence with 

For simplicity, let us assume for the moment that M is constant. Varying the action (2.12) we 
get the equation of motion for tt, 

4M\jf + SHtt) - {pq+pq)\t, - 3H{pQ +pq)t, = -\{pQ +pQ)h . (3.15) 

Small |c^| is equivalent to |pq +Pq| ^ so that the gradient and mass term in this equation can 
be neglected. Using again h = —26m [21], we thus have 

AM\if + SHtt) = {pQ + pQ)Sm . (3.16) 

One can verify that, neglecting decaying modes, the solution of this equation is 

^ . ,3.17) 

It is easy to show that this equation holds, for constant wq, for a general time dependent speed of 
sound which satisfies |c^(t)| <C 1. 

Equation (3.17) describes quintessence perturbations both for positive and negative 1 + wq. 
When WQ > — 1 quintessence energy density clusters in the dark matter potential wells, while in 
the opposite case wq < — 1 it escapes from them [27]. However, clustering of quintessence remains 
small compared to dark matter as the coupling with gravity is suppressed by 1 + wq . For very small 
values |1 + wq\ ~ M^/Mp^ we smoothly enter in the Ghost Condensate regime. Indeed it is easy to 
see that eq. (3.17) smoothly matches eq. (3.14) in the intermediate regime. 

It is important to stress that for wq < — 1 the speed of sound is constrained to be so small that 
quintessence effectively clusters on all scales. It would be interesting to understand the effect of the 
short scale clustering on structure formation. We will come back to this point in the conclusion, 
section 6. 

In this section we have studied the phenomenology of quintessence in various regimes of 1 + wq . 
Quintessence perturbations smoothly turn off when we approach the cosmological constant limit 
WQ = — 1 from both sides. This suggests that in general there is no pathology in crossing the 
WQ = — 1 line, as we discuss in the next section. 



4 Crossing the phantom divide 

It has been claimed that, during its evolution, single field quintessence cannot cross the wq = — 1 line 
as perturbations become pathological. For this reason this line has been dubbed "phantom divide" 
[7]. However, there is no real pathology in crossing this line, besides the fact that for wq < — 1 
short-scale gradient instabilities must be stabilized [8, 9]. If one does not take into account higher 
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derivative terms, a negative leads to catastrophic instabilities at short scales. Once instabilities 
are cured as we discussed in the previous sections, crossing the phantom divide becomes trivial.^ 

Indeed, the Lagrangian (2.13) gives an explicit way to construct a model which crosses the 
phantom divide. If one assumes for simplicity that quintessence is the only component in the 
Universe,^ the crossing of the phantom divide corresponds to a change of sign of H. In particular, 
considering only quintessence, one can use Priedmann equations to recast the Lagrangian (2.13) in 
the form (higher derivative operators will be considered later) 

P{X, ^) = -3M^,H\(I)) - Ml.HmX + 1) + ^MH</^')(X - 1)' • (4.1) 

This is similar to what happens in inflation, where the inflaton is the only relevant component in 
the Universe [2]. 

As an example, we consider the case where H evolves linearly in time and changes sign from 
negative to positive, 

m = ^{rnt-l). (4.2) 
This implies that H{t) will be a parabola of the form 

as shown in figure 3 (left panel). Using the general expression (4.1), we deduce that the Lagrangian 



P(<^,X) = -3 



^p^(f0'-0)+Mpii^J +,,\mc^-l)[[dct>f -l] + ^-M\ct>)[{dc^f + 1]" (4.4) 

admits the background solution (f) = t and the cosmological evolution (4.2) and (4.3). Note that 
there are no theoretical limitations on the choice of the background evolution H{t). Indeed, we can 
cross the phantom divide as many times as we want. For example, choosing H{t) oc sm{mt) the 
cosmological evolution keeps oscillating up and down around wq = — 1! 

We can now study perturbations around a solution crossing wq = — 1 to show that no pathology 
arises. The evolution equation derived from the action (2.12) reads 

(PQ +PQ +4M^) if+^dt [a^{pQ+pQ + 4M^)] ti-?,H{pq+pq)ti-{pq+pq)^ = -\{pQ+PQ)h . 

(X a Zi 

(4.5) 

At the phantom divide, the last three terms of this equation vanish but the equation is clearly non- 
singular. In our approach, it is manifest that the speed of sound squared changes sign at wq = — 1. 
In the example above is given by 

' ix\l - mt) + 2M4 ' ^^-"^ 



^Models that cross the phantom divide have been found in f{R) theories of gravity. However, for wq < — 1 

these models are equivalent to scalar fields with negative kinetic energy, i.e. to ghosts. Thus, according to 

our stability constraints, they are forbidden. 

^The conclusions drawn in this section also hold when including dark matter. 



16 




Figure 3: Example of phantom divide crossing, as given by eqs. (4-2) and (4-3), where we have 
defined He = — fi'^ / {2mMpi) . Left figure: behavior of H and H; the crossing of wq = —1 takes 
place at t = m^^ when H = and H = He- Right figure: trajectory on the quintessential plane. 



and the trajectory of the crossing on the quintessential plane is shown in figure 3 (right panel). The 
stability in the wq < — 1 region requires that |c^| remains extremely small so that it is mandatory 
to have a hierarchy between ^ and M, fi <C M. As we discussed, this hierarchy is naturally realized. 

Another approach to study perturbations is using a fluid description, as A:-essence is equivalent 
to a perfect fluid. In synchronous gauge, denoting with a prime the derivative with respect to the 
conformal time d7] = dt/a and defining TC = aH, the fluid equations read, in Fourier space, (see for 
example [28]) 

<5q + 3W(cI - n,Q)SQ = -(1 + u,q) [k^ + mHcl - c'j] ^ - (1 + u.^)^ , (4.7) 



^ + W(l-3cJ)^ = ^.«, (4.8) 

where 6q is the divergence of the velocity field of quintessence, 9q = ik'^T\/{pQ + pq) [21], is 
6pq/5pq calculated in a velocity orthogonal gauge {T^ = 0) [29] and it corresponds to the speed of 
sound squared that can be deduced from the vr Lagrangian'', to be distinguished from the adiabatic 
speed of sound squared, = Pq/pq = wq — wq/{3H{1 + wq)). The absence of pathologies at 
Wq — > — 1 can also be shown in this formalism. Indeed, in the continuity equation the divergence of 



,2 



^The velocity orthogonal condition = is equivalent to the condition tt = 0. As </) is unperturbed in 
this gauge, the perturbations of pressure and energy density only come from fluctuations of X, i.e. Spg/Spq = 
Pq,x /pq,x = Px / {'^PxxX + Px)i which is the speed of sound which appears in the tt Lagrangian (2.6). 
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is compensated by the prefactor in front of the squared brackets, while the 1 + wq term at the 
denominator in the Euler equation is harmless as also vanishes for wq —1. Thus both 5q and 
Oq are continuous through the divide. This is not surprising as 6q is just the Laplacian of the scalar 
perturbation it, = aQq/k'^ . 

At this point the reader may be puzzled: in the previous sections we stressed that close to the 
WQ = — 1 line quintessence behaves like the Ghost Condensate on cosmological scales, while eq. (4.5) 
as well as the fluid eqs. (4.7) and (4.8) do not contain higher derivative terms. Let us see why these 
additional terms are irrelevant in realistic cases of crossing the phantom divide. With these new 
terms, the equation of motion for tt derived from the full action (3.1) is obviously still continuous so 
that also in this case the line wq = — 1 can be crossed smoothly. The operator proportional to 
dominates in the Ghost Condensate strip around the wq = — 1 line. However, this happens only in 
the extremely narrow range |1 + wq\ < M^/Mp^ < 10~^^. The equation of state parameter wq will 
stay in this range only for a time much smaller than H~^, unless its evolution is tremendously slow. 
Thus TT has no time to evolve in the Ghost Condensate regime, so that for all practical purposes one 
can totally neglect this strip around the wq = — 1 line on cosmological scales. 

We have seen that A;-essence can be described with the fluid equations (4.7) and (4.8). Even 
including higher derivative terms, quintessence remains a perfect fluid (see appendix B) but does 
not satisfy the fluid equations (4.7) and (4.8) as these assume a linear dispersion relation. However, 
as we discussed, higher derivative terms are phenomenologically irrelevant on cosmological scales, so 
that one can still use the fluid description above when comparing with observations. 

Prom a practical point of view we conclude that, when comparing with observations a dark 
energy model which crosses the phantom divide, it is consistent and theoretically motivated to set 

= 0. On the other hand, it is inconsistent to turn off perturbations as sometimes done in the 
literature. 

5 Additional higher derivative operators 

As we discussed, higher derivative operators become relevant when the speed of sound is very close 
to zero. This regime is particularity interesting when wq < — 1 so that in the following we will 
consider mostly the case pq + pq < 0. 

Theories with very small should be thought of as tiny deformations of the Ghost Condensate 
theory [20, 1, 2]: in this limit one recovers the shift symmetry (j) ^ (/) + X, so that a small deviation 
from the Ghost Condensate is technically natural. In the Ghost Condensate limit there is an addi- 
tional symmetry that one can impose, i.e. the parity symmetry (p ^ —(p. The background = t in 
Minkowski space breaks this parity symmetry and the time reversal symmetry to the composition of 
the two; the theory of perturbations is thus invariant under vr —tt, t ^ —t [4]. This symmetry is 
present only when the background metric is Minkowski: in de Sitter there is a preferred time direc- 
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tion singled out by the expansion. In this case terms violating the symmetry will be proportional to 
H, and thus tipically suppressed by H/M. In this paper we have considered small departures from 
the Ghost Condensate limit, i.e., tiny breakings of the shift symmetry. These also generate terms 
which are not invariant under parity tt — — tt, t — > —t, as for instance the mass term in eq. (2.12). 
These terms will be of the same order of magnitude, i.e. suppressed by H/M as we are assuming 
that is the typical time scale of evolution of the operators. 

However, one can also consider the case when the parity symmetry — — is absent in the 
Ghost Condensate limit [1]. This happens for example if we add to the A;-essence Lagrangian (2.13) 
the operator 



C^, = - — (n(f> + 3H)iX-l), (5.1) 



M3 

Y 

which again docs not change the background evolution as it starts quadratic in the perturbations. 
For simplicity we assume that M = and a constant M.® In synchronous gauge at quadratic order 
this operator is 

£^ = M^TT + SHtt - SHtt - ^ + . (5.2) 

The first two terms in the parentheses contribute (after an integration by parts) to the time kinetic 
term. Assuming M ~ M they can be neglected in comparison with 2M^7r^. The third term gives a 
mass term that is parametrically smaller than H and can thus be neglected. 

To discuss the stability and phenomenology of this model, let us write the full action for pertur- 
bations, assuming \pq + Pq\ M^, 



-I 



r4 • 2 



1/ ^(Vv^)2 1, ,-V3. / h V^TT 



(5.3) 



where we have neglected the mass terms. This equation is analogous to eq. (3.1) for the M operator 
(2.16). Analogously to what we have done in sections 2.2 and 3.2 for the operator proportional to 
M^, we will now study the stability and phenomenology on cosmological scales with the operator Cj^. 
In appendix D we briefly study the effect of this operator at short distances, i.e. the modification 
of gravity induced by it. 

5.1 Stability constraints with M 

Let us first study the stability of the system neglecting other sources of gravity. The equation of 
motion for tt derived varying (5.3) reads 

^ ^ 3^ . _ PQ±|Q - ^Zjl = _PR±^k -^Ah + 3Hh) . (5.4) 



^This assumption can be relaxed by having a time dependence with time scale of order H ^, as in the case 

of 
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Notice that the operator M^'kV'^iT/a^ induces a spatial kinetic term for tt proportional to H. Indeed, 
this operator is a total derivative in Minkowski spacetime. Choosing M > 0, the spatial kinetic term 
has the "healthy" sign and can be chosen sufficiently large to cure the gradient instability associated 
to PQ +PQ < 0, giving a positive and very small c^. This also allows us to neglect the first term 
on the right hand side in eq. (5.4). To complete the stability analysis one has to take into account 
the mixing with gravity, i.e., solve for h in terms of the quintessence stress-energy tensor using the 
Einstein equation (2.23), and plug the result back in the right hand side of eq. (5.4). This, similarly 
to what happens for the Ghost Condensate, will give rise to a Jeans-like instability. 

The contribution to the stress-energy tensor of the operator (5.1) is (see appendix B) 

^PQ ^ (I - ' (5-5) 

6pQ D -2M^{TT + 2,H'k) . (5.6) 

Given the small speed of sound, time derivatives are much smaller than the spatial ones and the 
pressure perturbation is negligible, 5pQ <C Spq. Concentrating on frequencies much larger than the 
Hubble rate one can neglect the terms containing H-k and Hh in eq. (5.4). A further simplification 
comes from disregarding the standard A;-essence contribution to the energy density perturbation 5pQ^ 
i.e. 4M^7r, in comparison with M^V^7r/a^. Indeed, from eq. (5.4) we have M^V'^ir/a'^ ~ M'^ir/H » 
M^TT. Moreover, as we will see, the absence of Jeans instability will impose < Mp^H. This 
implies that the term with h in eq. (5.5) is negligible with respect to h in the Einstein equation 
(2.23), that becomes 

h = — ^ ^. (5.7) 

Plugging this into the right hand side of eq. (5.4) we finally find 

.. _ ( PQ+PQ , M^F _ \ VV _ 

This same result would have been found using a more rigorous approach, as done in [1]. Again, as 
in the Ghost Condensate case, mixing with gravity induces a Jeans-like instability, represented by 
the last term in this equation. Thus, for pq + pq < we need to cure both the gradient and the 
Jeans instabilities. This is possible for 

-{l + njq)^lq<^<l. (5.9) 

This stability window [1] is analogous to the one discussed in the Ghost Condensate case, eq. (2.29). 

We conclude that with the inclusion of the operator we can have a dispersion relation lu oc k 
with positive speed of sound squared; thus, there is no sign of instability even for pq + pq < 0. 
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5.2 Including dark matter 

Analogously to what done in section 3.2, to study the phenomenology induced by the M operator 
we study quintessence perturbations generated by the coupling with dark matter. For simplicity we 
assume matter dominance. The action (5.3) gives the equation of motion for tt sourced by the dark 
matter perturbation 

where we have used h = —2H5m and neglected the gradient term proportional to pQ + pq that is 
subdominant. Since the speed of sound = HM^/{4M'^) is very small and we are interested in 
cosmological scales, one would naively neglect the term with V^tt. This gives the solution 

However, the approximation of neglecting the gradients is not good. Indeed, when we plug this 
expression into the energy density perturbation 

6pQ=4M\ + M'(^^-^^ , (5.12) 

there is a cancellation of terms proportional to up to gradient terms. Thus one is forced to go 
back to the equation of motion (5.10) and keep the term proportional to the speed of sound. 
Once we do that, we obtain 

^\S^ + 1±^S„. (5.13) 



^ 24M4M2jnQ i?2a2 1-3u;q 

This equation displays the existence of two different regimes, in strict analogy with what happens 
in the Ghost Condensate case. For large enough |1 + wqI, the second term on the right hand 
side dominates and one recovers eq. (3.17), where the system behaves as standard A;-essence. Note 
however that even for pQ +PQ < there are no stability problems in the stability window (5.9). The 
dynamics of the system is dominated by the M operator only when we are very close to wq = — 1, 
i.e., for 

This region is the analog of the Ghost Condensate strip around the horizontal axis of figure 2. In 
this range the first term in eq. (5.13) dominates and Sq remains extremely small with respect to dm- 
Although the dispersion relation is of the form a; oc Cg/c, quintessence does not follow the simple fluid 
equations (4.7) and (4.8) because of the presence of higher derivative operators. However, as in the 
Ghost Condensate case, given the narrowness of the strip (5.14), for all practical purposes we can 
always use the fluid equations with = 0, even when crossing the phantom divide. In conclusion. 
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the addition of the operator M can stabihze A;-essence in the phantom case wq < —1, and the 
phenomenology of the model is the same as for a A;-essence with = 0. This general conclusion will 
hold even when considering both operators M and M at the same time, and can be extended to all 
the possible higher derivative operators, included in the general action in appendix C. 

6 Conclusion and outlook 

In this paper we have studied the most generic action describing the perturbations of a single field 
dark energy - here called quintessence - around a given background. We have constructed the 
action by adding to the A;-essence Lagrangian higher derivative operators that leave the background 
evolution invariant. Using this action, we have reproduced the results of [1] concerning the theoretical 
constraints on the equation of state parameter wq as a function of the speed of sound squared c^, by 
the requirement that perturbations are ghost-free - i.e., that their kinetic energy is positive - and 
that there are no gradient-like instabilities. These constraints have been conveniently represented 
on the quintessential plane (1 -|- wq)CIq vs. c^, in figures 1 and 2. 

In particular we have considered the case wq < — 1, which is commonly believed to be unstable, 
and we have shown that for very small both the gradient and the Jeans instabilities can be avoided 
and perturbations stabilized [1]. Higher derivative operators are crucial for the stabilization. Indeed, 
it is important to stress that taking an extremely small does not represent a fine tuning, as in 
the limit — we recover the Ghost Condensate theory which is protected by the shift symmetry 
(f) <p + \. Thus, for WQ < — 1 quintessence should be thought of as a small deformation of 
the Ghost Condensate limit [20, 1, 2]. When the higher order terms containing k'^ dominate over 
the spatial kinetic term c^k"^ the phenomenology reduces to that of the Ghost Condensate. This 
always happens on small scales, where the higher order gradients must dominate to stabilize the 
perturbations, but on cosmological scales this only occurs for values of wq extremely close to the 
one of the cosmological constant, i.e., for |1 -|- wq\Q,q < 10~^^. Away from this tiny strip - i.e., for 
all practical purposes - the behavior on cosmological scales is very different from that of the Ghost 
Condensate: higher derivative terms are irrelevant so that the phenomenology of the wq < — 1 side 
of quintessence reduces to that of a fc-essencc fluid with = 0. 

Furthermore, we have studied the behavior of quintessence perturbations when crossing the so- 
called phantom divide wq = —1. By restricting the analysis to fc-essence perturbations around a 
given background crossing the phantom divide, we have shown that, as the speed of sound vanishes 
exactly at the divide, perturbations remain finite during the crossing. For wq < — 1, higher derivative 
terms, while irrelevant on cosmological scales, are essential to stabilize the short scale gradient 
instabilities. We conclude that no pathology arises during the crossing: the phantom divide can be 
crossed without the addition of new degrees of freedom. We have illustrated this with an example 
shown in figure 3. An important thing to retain is that a consistent and theoretically motivated way 
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of comparing with data a dark energy evolution which crosses the phantom divide is to set to zero 
the speed of sound of perturbations. 

Our study motivates the possibihty that quintessence has a virtually vanishing speed of sound, 
especially when wq < — 1. Such a quintessence can be detected through its effects on structure for- 
mation. The speed of sound defines the sound horizon ig = a J Cgdt/a, which sets the characteristic 
length scale of smoothness of the perturbations. In the matter dominated era £q = 2csHq^ for a 
constant Cg. Hence, for Cg = 1 - corresponding to the speed of sound of a scalar field with a canonical 
kinetic term - quintessence can cluster only on scales larger than the Hubble radius while for = 
it clusters on all scales, thus affecting the gravitational potential and the formation of structures of 
dark matter and galaxies. The effect of a clustering quintessence can be measured with the cosmic 
microwave background [26, 27, 28, 30], galaxy redshift surveys [31], large neutral hydrogen surveys 
[32], or by cross-correlating the integrated Sachs- Wolfe effect in the cosmic microwave background 
with large scale structures [33, 34]. For instance, in [33, 31] it was found that with future surveys it 
will be possible to measure a zero speed of sound of dark energy if wq > —0.95. As these analysis 
were restricted to positive values of 1 + wq only, it would be interesting to repeat them for negative 
values. Notice that for a vanishing speed of sound, dark energy will actively participate to the for- 
mation of non-linear objects, affecting the halo bias. It would be interesting to evaluate this effect, 
which gives additional signatures of these models. 

Quintessence is perturbed by the presence of sources and thus modifies gravity as any other kind 
of matter. When this modification happens on scales much smaller than the Hubble radius one can 
properly talk about a theory of infrared modification of gravity. This happens when quintessence 
is close to the Ghost Condensate limit; in this case the modification of gravity is due to the Jeans 
instability induced by higher derivative operators, and persists even in Minkowski spacetime. Mod- 
ifications of gravity induced by the Ghost Condensate have been studied in details in [4, 22]. In this 
paper we have considered also the additional operator proportional to [1]; it would be interest- 
ing to investigate the deviation from General Relativity induced by this operator and its possible 
observational consequences, extending the preliminary analysis of appendix D, where the treatment 
has been restricted to linear perturbations in a Minkowski spacetime. 
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Appendix 



A Higher derivative operators in effective field theories 

In this appendix we want to study the Ghost Condensate and its deformations from an effective field 
theory point of view. In particular we want to show that, although the operator (V^tt)^ dominates 
the dynamics, operators containing higher time derivatives such as tt^ must be treated perturbatively. 
In the Ghost Condensate limit, the free tt action is 



S = — I df'xdt 



TT 



M2 



(A.l) 



where we neglected the mixing with gravity - as we are interested in the high energy behavior of 
the theory - and for simplicity we assumed that there is a single scale M (M 2± M). This action is 
manifestly invariant under the energy scaling [4] 

E^sE, t^s-H, x^s'^/^x, Tr^s^/\. (A.2) 

As the theory is not Lorentz invariant, time and space behave differently under rescaling, and tt does 
not scale as as in a Lorentz invariant theory. (See for instance [35] for an introduction to scaling 
in non-Lorentz invariant field theories.) 

What is the physical meaning of this scaling transformation? Assuming that the theory is 
weakly coupled, the free action gives the leading contribution to the correlation functions, so that 
these will be invariant under the scaling above. For instance, a relativistic massless scalar has scaling 
dimension 1. Thus, the two-point function satisfies^ 

(#)(x-y) = s-2(#)f^') ^ (00) oc^—. (A.3) 

\ s / |x — 

In the case of the action (A.l) above, the scaling transformation (A.2) yields 

(7r7r)(At, Ax) = s-'/'{nn) (^^, ^) . (A.4) 

Not only the scaling transformation gives information on the free theory, but, more importantly, 
it allows one to estimate the effect of different operators added to the free action. In particular, in 
the Ghost Condensate case, one can check that all additional operators allowed by symmetries have 
positive scaling dimensions, so that their importance is suppressed by E/M elevated to a positive 
power [4]. This implies that at low energy the theory is perturbative. For instance, the leading 
irrelevant operator is 7r(V7r)^, which has scaling dimension 1/4. 



^Note that the scaling dimension has nothing to do with the mass dimension of the field (p. Indeed, 
eq. (A.3) remains the same if we choose a non-conventional normalization of the action such that (f) has not 
mass dimension 1. 
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Operators containing higher time derivatives have positive scaHng dimensions so that they must 
be treated perturbatively. For instance vf^ has scaling dimension 2, so that at low energy it is 
negligible. Additional time derivatives naively suggest the existence of more and more solutions of 
the equations of motion. However, these solutions are non-perturbative in the expansion parameter 
E/M, and there is no reason to expect that they have any physical meaning. For example, taking 
seriously these solutions would imply that the Minkowski vacuum is unstable when considering 
higher order corrections to the Einstein-Hilbert action [36]. The correct way of treating these terms 
is perturbatively, i.e., evaluating them using the lower order equations of motion [37]. Following 
this logic, the additional solutions studied in the context of the Ghost Condensate theory in [38] are 
non-physical, as already pointed out in [39]. 

As we discussed in this paper, in certain regimes quintessence behaves as a deformation of the 
Ghost Condensate theory. The free action (A.l) is deformed by the addition of a (Vvr)^ term, 



S 



d xdt 



7^2 



M2 



(A.5) 



In these cases the dispersion relation is not exactly uj ~ I? jM but it contains also a linear term 
io ~ Cs/c, with Cg <^ 1. (For this discussion we assume that is positive.) 

The situation is now trickier than before because one cannot find a scaling transformation which 
leaves the full action (A.5) invariant. On the other hand, one can separate two regimes, depending 
on which of the gradient terms dominates, as illustrated in figure 4. 



k 

cutoff 



c,M 







GC: u - k'^/M 



uo ^ c.k 



Figure 4: The two scaling regimes as a function of the momentum k together with the scaling 
dimensions of some of the operators. Above: the Ghost Condensate regime where uj ~ k'^/M; below: 
the regime where uj ~ Csk. 



For k ^ CgM, the dispersion relation is dominated by (V^Tr)^ and the theory behaves as the Ghost 
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Condensate. In this regime the scaling of all additional operators can be obtained from eq. (A.2). 
Notice that now there is a relevant operator, (Vtt)^, that becomes more and more important at 
low momenta (and energies). The coefficient of this operator is however suppressed by the small 
deformation parameter Cg ^ 1. Thus, it can be treated perturbatively as long as A; S> CgM. 

On the other hand, when k <^ CgM the (Vtt)^ operator dominates the free action. In this regime 
the scaling becomes the same as in the relativistic case. Now both tt^ and (V^tt)^ are irrelevant 
operators with the same scaling dimension 2. However, time and spatial derivatives are still on a 
different footing because the time derivatives are suppressed by Cg with respect to the spatial ones, 
u ~ Cgk. Thus TT^ ~ Cg{V^TT)^ for k ^ CgM. In the intermediate regime k ~ CgM this suppression 
can also be obtained from the Ghost Condensate limit. Indeed, at high momenta the operator tt^ 
scales like (k/M)^, so that it is suppressed by Cg for k ~ CgM. 

Even though these theories make perfect sense as effective field theories, it is not clear whether 
one can find a UV completion. In particular, the violation of the null energy condition may be 
problematic in the context of black hole thermodynamics [40] . 

B Stress-energy tensor and fluid quantities 

Here we compute the stress-energy tensor for the fc-essence action with the addition of the two higher 
derivative operators Cj^^ and Cj^, i.e.. 



4 / 



P{^, X) - ^{D<l> + 3H{^)f - ^{D<l> + 3H{ct^)){X - 1) 



(B.l) 



We start by computing, using eq. (2.7) the stress-energy contribution of the fc-essence action, which 
can be written as 

(T^,)p = 2Px{(t>,X)d^<t>du4> + Pi^,X)g^, . (B.2) 
It is well known that this stress-energy tensor can be put in the perfect fluid form, 

Tf,u = (p + p)Uf,Uu + P9f,u , (B.3) 
by defining the /c-essence contribution to the energy density and pressure respectively as 

pp = 2XPx-P, (B.4) 
PP = P, (B.5) 

and the unit four- velocity of an observer comoving with the fluid as 

u, ^ -M= ■ (B-6) 
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The stress-energy tensor (B.2) can be expanded in perturbations around the background (f){t) = t. 
In synchronous gauge one finds 



(Too)p = PQ + PQT^ + {pQ +PQ + 4M^)7r , 
{Toi)p = {pQ+PQ)diTr, 

(Tij)p = a^bij \pQ+PQn + {pQ +PQ)n] + a^pqhij . 



(B.7) 
(B.8) 
(B.9) 



One can do the same for the operators proportional to and and show that they can also 
be written in the perfect fluid form (B.3), by noting that defined in eq. (B.6) can be written in 
several different ways, 

d^U(t> _ d^{U(t> + ?,H{ct))) 



^^{dxf ^-{du^f ^-[d{u^ + m)f 



(B.IO) 



Let us first compute the contribution to the stress-energy tensor of the operator proportional to M^, 
which reads^° 



1 



29(^(n</.+3i?)5,) - (□</.-3i?) (□</.+3ff)5^, (□</)+3i/)a/3</)5, 



, (B.ll) 



where (^j,) denotes symmetrization in fi and v. The perfect fiuid form (B.3) can be obtained by 
recognizing its contribution to the energy density and pressure, defined as 



Pm2 = 



Pm2 = 



g'^'d^iDcf) + 3H)d^(l> --(□</)- SH){n(f) + 3H) 
g^'^df.incp + 3H)d^(P + ^(00 - 3H)in4> + 3H) 



In order to expand in perturbations let us first compute the Laplace-Beltrami operator, 
n4> = -La^ {^gt'''d,4>) = -3H -jr-3H7r + ^-h. 



-9 



Then the stress-energy tensor becomes 



(B.12) 
(B.13) 

(B.14) 



{dt + 3H) + 3Hii - 3Hn 



V^vr 1 



V^TT 1 ; 



+ 25o(z.5^) 7f + 3H7r - 3Ht: ^ + -h 



(B.15) 



We now perform the same procedure for the operator proportional to . Its stress-energy 
tensor is 



-(□</> + 3H)d^<pd,(P - d^^Xd,)(t^ - ^9^u{X - 1) + ^5M^a"Xa„<; 



(B.16) 
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Here we take for simplicity M and M to be constant. 
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and once again it can be written in the perfect fluid form by defining 



Pm3 
Pm3 = 



^g^'^'d^xd^ct^ - {u<t> + m)x + Ih{x - 1) 



2 

^-g^-d^XdA-^-H{X-l) 



(B.17) 
(B.18) 



When expanded to first order in the perturbations, the stress-energy tensor for the operator propor- 
tional to takes the form 



OOjMS 



M 



3 ( eifTT - m-K 



3 



V^TT 1 



(B.19) 

(B.20) 
(B.21) 



Using the notation of [21], the stress-energy tensor of a perfect fiuid in synchronous gauge can 
be written as 



(B.22) 
(B.23) 
(B.24) 



Too = PQ + 
Toi = -a{pQ+pQ)d^V~^eQ, 



Notice that we defined 9q in the main text, below eq. (4.8), using conformal time [21]. Here we 
are in cosmic time and this introduces the factor of a in the equation for Toi above. Thus, from 
equations (B.7)-(B.9), (B.15), and (B.19)-(B.21) one gets, for our complete stress-energy tensor. 



^PQ =PQ'^ + {PQ +PQ + ^M'^)-k 



+ {dt~3H)(^jf + 3H7r-3H7T- 

-^.f ■ V^TT l-\ 



V^TT 1 ■ 

^ + 2^ 



V^7r + 



PQ + PQ 
SPQ =PQT^ + [PQ +Pq)'^ 



7f - SHtt + SHtt + ^ --h] + V^^ , 

PQ + PQ 



(B.25) 
(B.26) 



+ M2(5t + 3H) + 3Hi: - SHir - ^ + ^hj - M^{dt + 3H)tt . (B.27) 



C General action in unitary gauge 

We wish to write the general action for the perturbations of a single quintessence field minimally 
coupled to gravity with no direct couplings to other fields. Following ref. [1, 2], we choose a gauge 
where the scalar field perturbation is set to zero but it appears as a scalar metric degree of freedom. 
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In this gauge the constant time hyper surf aces are equivalent to the uniform field hyper surfaces. In 
the fluid language, this implies that the velocity is orthogonal to the constant time surfaces, = 0, 
that is why this gauge is called 'velocity orthogonal' [29]. Using a particle physics terminology we 
can also call it 'unitary gauge', as all the degrees of freedom are in the metric. Notice that one can 
always parameterize the perturbations of the scalar field as 

cj)it,x) = Mt + <t,x)), (C.l) 

where tt is the difference between the constant time and uniform field hyper surf aces, so that the 
unitary gauge corresponds to tt = 0. 

In order to find the effective action for quintessence, gravity and other matter components 
described by the Lagrangian Cm, we write down all the terms that preserve the symmetries of the 
system. Our choice of gauge breaks time diffeomorphism invariance while preserving invariance under 
spatial diffeomorphisms. Thus, we include linear and quadratic combinations of generic functions of 
time t, the time-time component of the inverse metric and the extrinsic curvature of the constant 
time hyper surf aces. The effective action up to second order in perturbations is [1, 2]: 



/ 



d^x dt 



(0.2) 



where R is the Ricci scalar and Kij is the extrinsic curvature of constant t hypersurfaces, which at 
linear order reads 

1 



Kij = - a/ -5°° [dQQij - digoj - djgio) , (C.3) 

and we have defined 6Kij = K^j — a^HSij, and SK = K\ — 3H. 

There are other operators that are invariant under spatial diffeomorphisms that one would naively 
include in this action. However, these operators are irrelevant at energy scales below the cutoff 
M ~ M ~ M. For instance, one could include the operator (5™)^. This would give the term tt^ 
in the final action, which indeed also appears as part of the operator (2.16), once expanded in the 
perturbations, eq. (2.17). However, as explained in section 2.2, for frequencies smaller than the 
cutoff TT^ is negligible with respect to (V^tt)^ so that it can be ignored in the action. For simplicity 
we will also ignore the operator proportional to as it leads to terms qualitatively similar to those 
proportional to M^. 

One can easily fix the coefficients c{t) and A(t) by computing the background stress-energy 
tensor. This gives 

PQ = A{t) - c{t) , (C.4) 
PQ = -c{t)-K{t). (C.5) 
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Finally, using these relations, we can write the action (C.2) in unitary gauge in terms of the back- 
ground quantities PQ{t) and PQ{t), 



-I 



d X^/^ 



^R + Cm +PQ - 1{PQ+PQ){9'' + 1) + ^(5°° + 1)^ 



SK{g^ + 1) 



(C.6) 



Now we want to rewrite this action in a gauge-invariant form. This can be done by performing 
the following time-coordinate transformation, 



X X = X 



t^t = t + 7:{x) 

that reintroduces tt defined in eq. (C.l). The action for tt reads, up to second order. 



(C.7) 



^= J j PQ + PQTT -I- ^PQTT^ 



(5°° + 1) - 2^ + 2(5°° + 1)71 - TT^ + 2g'%7r + ^ 



(C.8) 



[(^ + 1) - 27r; 



5K - SHtt - [(5°° + 1) - 27r' 



while the part of the action containing R and Cm is invariant under general diffeomorphisms. We 
now choose to work in the synchronous gauge, which is defined by the metric 



ds^ = -dt^ + a^{t){5ij + hij)dx'dx^ 



(C.9) 



Using the notation of [21], the two scalar degrees of freedom of hij are its trace h = S^^hij and rj 
which is defined by V^/iy = didjh + 6(didj — ^6ijV^)r]. Using this metric, after integrating by parts 
and using the background continuity equation pQ + SH{pQ + pq) = 0, the action (C.8) takes the 
form 



1 



S = I d xa \PQ + ^(pq+Pq) 



(Vtt) 



+ 2MV + ^H{pQ+pQ)Tr^ - ^{pQ+PQ)hTr 



(C.IO) 

We can now compute the stress-energy tensor of quintessence using the action (C.8). Expanding 
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in the perturbations, its components read 



Too = PQ + [PQ +PQ + 4M^) + PQTT 

- 3HM^ (\h - SHtt - ^TT^ + M" 
Toi = {pQ + PQ)diTr + M'^di ( h - 3H7T - 



2 

V^TT 



Tij = pga^Sij + \pQ7T + {pQ + pq) tt] a^dij + pqa^hi 



1- . V2 

h — 3/f7r ^TT + Si^TT 



+ 2M^a^%(5o + 3H) -h - SHn - 



-2M^a^Sij{do + 3H)'k. 



(C.ll) 
(C.12) 

(C.13) 



In the main body of this paper we have constructed the action for vr using a different procedure 
from the one presented here. Indeed, we started from the action of fc-cssence, eq. (2.6), and we 
added the two (/^-dependent higher derivative operators — M'^ + 3H {(p)]'^ /2 and —M'^((p)[D(j) + 
3H{(I))]{X — l)/2, that do not change the background equations of motion. Also the action (C.IO) 
can be constructed similarly. First of all, note that the first line of eq. (C.IO) is the action for 
fe-essence. Indeed, it is equivalent to the action (2.12), which was found by expanding the fe-essence 
action (2.1) in terms of tt. The second line of eq. (C.IO) can be constructed by noting that the 
extrinsic curvature of the hypersurfaces of constant 4>, defined as 



where is the unit vector orthogonal to (p defined in (B.6), can be rewritten as 

1 



(C.14) 



-{d<i>) 



2 V p^j/v 



(C.15) 



When expanded around the background solution (/) = t, then SK = K'^n — 3H{(f)) reads^^ 



V2 1 . 
5K = -3H'K ^7r + -/i. 



(C.16) 



Thus the action (C.IO) can be constructed by simply adding to the fc-essence action the two operators 
-M'^{(f))5K'^/2 and -M^{4>)5K{g^^ + l)/2, that do not change the background solution. 



D Modification of gravity with M 

In the main text we studied the effects of the operator M, focusing on stability and on the phe- 
nomenology at cosmological scales. In analogy to what happens for the Ghost Condensate, we expect 

^^One may wonder why terms of the form ttV^tt do not appear in 5K^, while they do appear in {D(f> + 
3H)^. As seen in cq. (C.15), a time diffeomorphism docs not change the extrinsic curvature of constant (p 
hypersurfaces. Thus, SK does not contain n and ttV^tt is not generated by SK"^. 
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that this operator will also be relevant at short scales, inducing a modification of gravity. In this 
appendix we perform a preliminary analysis, restricted to linear perturbations only, although the 
non-linear dynamics has been shown to be relevant and quite rich in the Ghost Condensate case 
(see for example [22]). To simplify the analysis we set pQ + PQ = and M = 0. Although the 
background quintessence stress-energy tensor is the one of the cosmological constant, there is still 
a propagating scalar degree of freedom. Its mixing with gravity induces a deviation from General 
Relativity; indeed the Ghost Condensate was originally proposed as a consistent modification of 
gravity in the infrared. The simplest setting to study this modification of gravity is in the Newto- 
nian regime u'^ <^ around Minkowski spacetime, where the new scalar degree of freedom modifies 
the Newtonian potential For this purpose we will closely follow the discussion done in [4] for the 
Ghost Condensate case, i.e., for the operator M. 

Working in Newtonian gauge with ^ = ^, the metric is ds^ = — (1 + 2$)dt^ + (1 ~ 2^)dx^ and 
the quadratic Lagrangian for tt and $ reads 

C = -M|i(V$)2 + 2M^($ - irf + m3($ - ^)(4$ - vf + VV) . (D.l) 

Let us first assume that M is time independent. Dropping total derivatives and terms which are 
negligible in the limit uP <^k'^, we are left with 

C = -M|i(V$)^ + 2M'^($ - -kf + mVV^$ . (D.2) 

In terms of the canonically normalized fields tTc = 2M^7r and $c = \/2Mpi$, the Lagrangian in 
Fourier space can be written as 



with 



Setting to zero the determinant of this matrix gives the dispersion relation 

which reproduces the Jeans instability already shown in eq. (5.8). The Jeans instability arises from 
the non-diagonal (mixing) term and it is thus proportional to instead of M^. 

To study the corrections to the Newtonian theory, one can look at the propagator of $ that is 
the (4>, entry of A^~^. This can be written as 



1 



k'^M^ 1 



(D.6) 
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where the term — is simply the standard Newtonian propagator. As expected, a substantial 
deviation requires, for a given distance a sufficient time for the Jeans instability to develop, 

- *° k\ (D.7) 



We can now consider the case of a time dependent M. In this way we introduce new terms in 
the action that were previously dropped because total derivatives. The same happens if we had 
considered a time dependent spatial metric, but here we stick to Minkowski for simplicity. The 
matrix M becomes 

where H is the typical rate of variation of M^, = HM'^ (if the time dependence is induced by 
the metric this becomes the Hubble rate). Computing the determinant and restricting to frequencies 
much larger than H we get the dispersion relation 

which correctly matches eq. (5.8). The value of M can be chosen to avoid the Jeans instability and 
have a healthy dispersion relation. The propagator becomes 



1 



^ k'^M^ 



(D.IO) 



The scalar degree of freedom induces a 1/r force which adds to the Newton law: this force, 
however, propagates at a very small speed 

Given the absence of a Jeans instability, the modification of gravity induced by M is very different at 
linear and non- linear level with respect to the Ghost Condensate case [4, 22]. More work is needed 
to understand the constraints on M coming from the modifications of gravity that it produces. 
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